Universality in low energy states: from few-atoms to few-nucleons 
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We investigate universal behavior in elastic atom-dimer scattering below the dimer breakup thresh- 
old calculating the atom-dimer effective-range function ak cot 5. We compare our results with the 
universal zero-range form deduced by Efimov, ak cot S — ci{ka) + C2(fca) cot[so ln(aK,) -I- (j}{ka)], 
for different values of the two-body scattering length a and of the three-body parameter k,. We 
observe a good agreement introducing a particular type of finite-range corrections. Furthermore, 
we show that the same parameterization describes a very different system: nucleon-deuteron scat- 
tering below the deuteron breakup threshold. Our analysis reveals a universal behavior that ranges 
from few-atom systems to few-nucleon systems, and clarifies the nature of the pole structure in the 
effective-range function of nucleon-deuteron scattering. 



C/2 



I 

o 



> 
m 



X 



Introduction. Scattering of two particles at very low 
energy shows universal behavior encoded in the scatter- 
ing length a and in the effective range Vg. In fact, sys- 
tems with different interactions sharing the same scat- 
tering length and the same effective range have the same 
effective range function, kcotS = — 1/a -I- and, 
accordingly, the same low energy behavior. In the limit 
a 3> ^0 , where the scattering length is much greater than 
the typical range tq of the potential, not only the scat- 
tering process is universal, but also some bound-state 
properties. In fact, when a +oo (known as unitary 
limit), the two-particle system has a shallow-bound state 
with the bound-state energy E2 ~ fi^/ma?' fixed by the 
scattering length. In this limit, the physics is scale in- 
variant. 

In the 1970s, V. Efimov d [3] showed that the scale 
invariance is broken in the s-wave three-body sector of 
a bosonic system. The residual symmetry is the discrete 
scale invariance (DSI); namely, the physics is invariant 
under the rescaling r — >■ A"r, where the constant is usu- 
ally written A = with sq ~ 1.00624 an univer- 
sal number that characterizes a system of three-identical 
bosons. One consequence is that at the unitary limit 
the three-body spectrum consists of an infinite number 
of states that accumulate to zero with the ratio between 
two consecutive states being E^~^^/E^ = q-'^'^/^o_ por 
finite scattering length, the binding energies satisfy the 
Efimov's equation 
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with tan^ = -{mEr^ /h^f/'^a. The function A(0 is uni- 
versal and a parameterization in the interval [— tt, ~7r/4] 
is given in Ref. [3|. The three-body parameter the 
wave number of the n = state at the unitary limit. 

The DSI constrains the form of the observables to be 
log-periodic functions of the control parameters. One 
example is the atom-dimer scattering length which has 



the general form 

a AD /a = di+ d2 tan[so ln(aK* ) + ^3] , (2) 

where di,d2, ds are universal constants to be determined. 
For collisions below the dimer breakup threshold, DSI im- 
poses the following universal form for the effective range 
function 

ka cot S = ci{ka) + C2{ka) cot[so ln(aK*) + (j^ika)] (3) 

with 5 the atom-dimer phase-shift and ci, C2, (/> universal 
functions of the dimensionless variable ka, where k"^ — 
{A/?>)E/{h'^ /m), being E the center of mass energy of 
the process. As fc — >■ 0, fcacot^ — ?► —a/aAD and at fc = 
the constants di,d2,dz and ci(0), 02(0), 0(0) are related 
by simple trigonometric relations. 

In this paper we study in detail the universal behav- 
ior of a AD and of the effective range function ka cot 5. 
To this aim we use the family of atomic ^He-^He poten- 
tials derived in Ref. 4] for several values of o, running 
from a « 440 a.u. to a w 50 a.u.. The corresponding 
dimer energies range from E2 ~ 0.22 mK to E2 ~ 21 
mK covering two order of magnitude. For selected val- 
ues of a in the mentioned interval we calculate a ad and 
the s-wave phase-shifts in order to construct the effective 
range function below the dimer breakup threshold. 

Our numerical results are used to analyze the univer- 
sal form of Eqs. ([2]) and ([3]). As the value of a increases 
we observe that a ad changes sign tending to —00. This 
behavior produces a pole in the effective range function. 
Furthermore, as the calculations are done using finite- 
range interactions, we can extract finite-range corrections 
by comparing our results to the zero-range theory. In- 
terestingly, for the explored zone of positive scattering 
length, the range corrections can be taken into account 
by a shift in the variable K^a. 

Finally, the universal character of the effective range 
function can be used to evaluate a very different sys- 
tem: low energy nucleon-deuteron scattering. It is well 



known that the nucleon-deuteron effective range func- 
tion presents a pole structure that has been related to 
the presence of a virtual state [5]. We shown that this 
structure is related to the universal form given by Eq.Q 
and, using the parameterization determined in the atomic 
three-helium system, we shown that this equation can be 
used to describe nucleon-deuteron scattering as well. In 
this way, the universal behavior imposed by the DSI is 
analyzed in systems with natural lengths that differ of 
several order of magnitude. 

The three-boson model. We construct the model using 
the LM2M2 [||, one of the most used ^He-^He poten- 
tials, as the reference interaction, with the mass param- 
eter h^/m = 43.281307 (a.u.)^K. Following Ref. fl we 
define an attractive two-body Gaussian (TBG) potential 

V{r) = Vo e-'-'/^o , (4) 

with range vq = 10 a.u. and strength Vq fixed to repro- 
duce the values of a given by V\ (r) = A • T4m2M2 (r) ■ The 
unitary limit appears at A ~ 0.9743 whereas for A = 1 
the values of the LM2M2 are recovered: a — 189.41 a.u., 
S2=-1.303 mK and re = 13.845 a.u.. 

The use of the TBG potential in the three-atom sys- 
tem produces a ground state binding energy appreciable 
deeper than the one calculated with V\{r). For example, 
at A = 1 the LM2M2 helium trimer ground state binding 
energy is = 126.4 mK whereas the one obtained using 
the two-body-soft-core potential in Eq. ^ is 151.32 mK. 
Hence, we introduce a repulsive hypercentral-three-body 
(H3B) interaction 

VF(pi23) = Wo e-''-3/Po , (5) 

with the strength Wq tuned to reproduce the trimer en- 
ergy i?3 obtained with Vx{r) for all the explored values 
of A. Here pfja = |(^12 + ^23 + ''31) hyperradius of 

three particles and pn gives the range of the three-body 
force. Following Ref Q we use po = ''o- 

Varying A from the unitary limit to A = 1.1 we obtain 
a set of values for the ground state binding energy 
and first excited state E} using the TBG and TBG+H3B 
potentials in a broad interval of a. We can use the results 
for i?3 at the unitary limit (using Eq.([T])) to determine 
= 0.002119 a.u. and = 0.001899 a.u. for the TBG 
and TBG+H3B, respectively. In Fig. [1] we collect our 
results for the ratio E\/E2 as a function of K*a for the 
TBG potential (circles) and of the TBG-hH3B potential 
(squares). We can compare the numerical predictions to 
the predictions of the Efimov's binding energy equations 
of Eq.([T]), given in the figure by the dashed line. We can 
observe that the numerical results lie on a curve shifted 
with respect to the dashed line. We can interpret the 
shift as a consequence of the finite-range character of 
the numerical results. Accordingly, we can adapt the 
Efimov's equation to treat finite-range interactions. Once 
we fix n = for the first excited state E\, we rewrite 
Eq. ([D) as follow 

Ell E2 ^ t&rv^ S, , <a = exp[-A(0/2so]/cose, (6) 




FIG. 1. (Color online) Energy of the first excited state of 
the trimer as a function of K^a. The dashed Une is the uni- 
versal prediction of the Efimov law, while the solid line is 
the translated universal curve. The circles and squares are 
the calculations using the TBG and TBG+H3B potentials 
respectively. 

where we have introduced a finite-range three-body pa- 
rameter 

k',-kJ1 + € + ■■■), (7) 

with e = r*/a. Therefore, and k'^ are equal at the uni- 
tary limit (a 00) or in a zero-range theory (r* — s> 0). 
In fact, the usual form of the Efimov's equations is recov- 
ered in both limits where E2 goes to zero (unitary limit) 
or E2 — ^ h^/rna^ (the zero-range limit). In the present 
analysis we find = 21 a.u. ~ 2rQ. The corresponding 
results, obtained by solving Eq. ([6]) with 
are shown in Fig. [T] as a solid line. 

atom-dimer scattering. To describe atom-dimer scat- 
tering below the dimer threshold we calculate aAD and 
the s-wave atom-dimer phases 6 using the TBG and 
TBG-I-H3B potentials at different energies. We use the 
hyperspherical harmonic (HH) method in conjunction 
with the Kohn variational principle Applications of 
the method to describe a three-helium system with soft- 
core interactions as used here can be found in Ref. [8j- 
Defining E2 = h^/ma^, in Fig. [5] we show the results for 
the ratio clad/ 0,3 in terms of the product K*a. It can be 
observed that the calculated points, given as full circles 
(TBG potential) and full squares (TBG-hH3B potential), 
lie on a curve shifted with respect to the dashed line rep- 
resenting Eq.® with the parameterization of Ref. [3]. 
We can interpret again the shift as produced by the finite- 
range character of the calculations. Accordingly we can 
adapt Eq.Q to describe finite-range interaction as 

o.ad/o.b = di + d2 tan[so In(K'^a) + ds] . (8) 

In fact, replacing in the above equation k'^u — K*a-f K*r*, 
with the same numerical value of before, the solid 
line is obtained in Fig. [21 The values di = 1.531 
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d2 = —2.141, da = 1.100 slightly modify the parameteri- 
zation of Ref. [3] to better describe the numerical results. 
This new parameterization is shown as a dotted line in 
the figure where we can observe an improvement in the 
description of the results close to the unitary limit. 

The shifted formula can be used to determine the ra- 
tio a*/a_, where a_ is the scattering length at which the 
three-body states disappear into the three-atom contin- 
uum, and a* is the scattering length at which the three- 
body states disappear into the atom-dimer continuum. 
For the potential models used in this work the values of 
a_ are given in Ref. Q whereas the values of a* can be ex- 
tracted by equating the argument of the tangent in Eq.® 
to —tt/2. Using these inputs we obtain a^/a^ w —0.4. 
The zero-range universal formulas Eq. ([IJ and Eq. ([2]) 
predict a*/a- = —1.06, but recent experimental results 
give lower values for this ratio [§-[ll| ; namely, it has been 
measured a*/a_ ~ —0.7 in potassium [oj and in bosonic 



lithium 1 1 



while a^/a- « —0.4 in cesium 11 ul . The 
difference is given by finite-range effects |12Hl6j . which 
in our case are encoded in the shift. 
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FIG. 3. (Color online) The ratio aAo/ciB as a function of 
EI/E2. Open circles and open squares correspond to TBG 
and TBG-I-H3B potentials respectively. The dashed line cor- 
responds to Eq. ([2} with the parameterization of Ref. 
whereas the dotted line shows the present parameterization 
of Eq. ©. 
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FIG. 2. (Color online) Universal plot for aAo/aB in terms 
of K*a. Open circles and open squares correspond to TBG 
and TBG+H3B potentials respectively. The dashed line cor- 
responds to Eq. ((21), whereas the solid line corresponds to 
Eq. ((8]). The dotted line shows the present parameterization 
of Eq. (O. 



It is interesting to see that the finite-range corrections 
cancel in the description of the scattering length as a 
function of the trimer energy. This is shown in Fig. [3] 
in which the present calculations and the zero-range uni- 
versal theory, Eqs. ([1]) and ([2]), are in close agreement. 

We now present results for atom-dimcr scattering at 
energies below the dimer breakup threshold for differ- 
ent values of a. We adapt Eq. ([3]) to finite-range in- 
teractions by considering kas cot S as a function of the 
dimensionless center of mass energy (fca)^ and using 
At'^a = K^a + Ki,r,, in the argument of the logarithm func- 



kas cot d — ci(fca) + C2{ka) cot[so ln(aK'^) + (pika)] . (9) 

In Fig. 13] we show our results (given as full squares) at 
different values of tt'^a. In the figure we can observe very 
different patterns. For the smallest values of K'^a the 
behavior is almost linear in all the energy range. Start- 
ing at values of n'^a « 0.4 a curvature appears close to 
zero energy, pointing out to an emergent pole structure 
that becomes evident at larger values of K^a. Specifically, 
the pole appears when qad changes sign (see Fig. [2]) or, 
as given in Eq.([9]), when the argument of the cotangent 
function becomes zero (or rnr). The shadow plot in the 
first row of Fig. [4] corresponds to the case A = 1 and de- 
scribe '*He-'*He2 scattering (full triangles). The shadow 
plot in the second row corresponds to nucleon-deuteron 
scattering as discussed below. The solid curves are ob- 
tained using the finite-range-adapted Eq. (|9]) with the pa- 
rameterization of Ref. We can observe a noticeable 
agreement along the whole range of values. 

nucleon-deuteron scattering. The universal effective 
range function has been determined using the TBG and 
the TBG+H3B potentials describing an atomic three- 
helium system. The universal character of the function 
allow us to apply it to describe a very different system: 
nucleon-deuteron (n — d) scattering. To this aim we use 
the results of Ref. Q in which n — d scattering has been 
described using a spin dependent central potential. In 
that reference they obtained and — 0.71 fm for the n — d 
scattering length, and the effective range function has 
been parameterized as 



k cot 6 
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1 + Ec.m./ Ep 



(10) 



with Eri 



T60 keV and 



T27 fm. The pole ap- 
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FIG. 4. (Color online) The effective range function at different values of K»a as a function of (ka)^. The square points are our 
calculations; the triangle points are the calculations for real '*He-'*He2 scattering. The full circles, for K,a = 0.543055, correspond 
to neutron-deuteron scattering in the doublet channel. The solid curves have been calculated using the translated-universal 
formula for the different values of «:»a. 



pears in the negative region similar to what happens in 
Fig. |4] at intermediate n'^a values. From the values of 
Ep and and it is possible to determine the corresponding 
values of a and k'^ by using the universal function (j) and 
Eq. dSl). We obtain a = 4.075 fm and <a = 0.5779. The 
shadow panel of the second row in Fig. |4] shows a com- 
parison of the universal function (solid line) to the n — d 
scattering results of Ref. "E] (full circles). We observe a 
noticeable agreement. It should be noticed that the spin 
dependent potential used in Ref [1| reproduces the singlet 



(^a„p « —20 fm) and triplet {'^Unp ~ 5 fm) n—p scatter- 
ing lengths. Accordingly, the three- nucleon systems has 
a symmetric plus a mixed symmetry component. The 
value of a extracted from the universal function, which is 
close to ^a„p, can be considered an effective value in an 
equivalent three-boson system with the given and and the 
corresponding effective range function. A more deeper 
analysis extending the model to spin dependent interac- 
tions is at present underway. 

In conclusion we have calculated the effective range 
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function in three-boson systems using the three-hehum 
system as a reference. We have shown that the results 
can be described by the zero-range universal formula us- 
ing a shifted three-body parameter n'^a = K^a + n^tr^,. 
Interestingly, the value of r, necessary to describe the 
results is the same for £'3, for a ad and for the effective 



range function. Furthermore, we used this function to 
describe low energy n — d scattering finding a quantita- 
tive agreement with calculations on this system. This 
analysis connects the universal behavior of atomic sys- 
tems having large two-body scattering length to nuclear 
systems. 
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